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ABSTRACT 

We present a multiwavelength morphological analysis of star forming clouds and filaments in 
the central (<: 50 kpc) regions of 16 low redshift (z < 0.3) cool core brightest cluster galaxies 
(BCGs). The sample spans decades-wide ranges of X-ray mass deposition and star forma¬ 
tion rates as well as active galactic nucleus (AGN) mechanical power, encompassing both 
high and low extremes of the supposed intracluster medium (ICM) cooling and AGN heating 
feedback cycle. New Hubble Space Telescope (HST) imaging of far ultraviolet continuum 
emission from young (<; 10 Myr), massive (;> 5 M 0 ) stars reveals filamentary and clumpy 
morphologies, which we quantify by means of structural indices. The FUV data are compared 
with X-ray, Lya, narrowband Ha, broadband optical/IR, and radio maps, providing a high 
spatial resolution atlas of star formation locales relative to the ambient hot 10 7-8 K) and 
warm ionised (~ 10 4 K) gas phases, as well as the old stellar population and radio-bright AGN 
outflows. Nearly half of the sample possesses kpc-scale filaments that, in projection, extend 
toward and around radio lobes and/or X-ray cavities. These filaments may have been uplifted 
by the propagating jet or buoyant X-ray bubble, or may have formed in situ by cloud collapse 
at the interface of a radio lobe or rapid cooling in a cavity’s compressed shell. Many other ex¬ 
tended filaments, however, show no such spatial correlation, and the dominant driver of their 
morphology remains unclear. We nevertheless show that the morphological diversity of nearly 
the entire FUV sample is reproduced by recent hydrodynamical simulations in which the AGN 
powers a self-regulating rain of thermally unstable star forming clouds that precipitate from 
the hot atmosphere. In this model, precipitation triggers where the cooling-to- freefall time 
ratio is t coo \/tff ~ 10. This condition is roughly met at the maxmial projected FUV radius for 
more than half of our sample, and clustering about this ratio is stronger for sources with higher 
star formation rates. 

Key words: galaxies: active - galaxies: star formation - galaxies: clusters: intracluster 
medium - galaxies: clusters: general 
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1 INTRODUCTION 

Many giant elliptical, groups, and clusters of galaxies inhabit an 
X-ray bright halo of ;> 10 7 K plasma whose core radiative life¬ 
time is much shorter than its age. Absent a heating mechanism, 
simple models predict that the rapid cooling of gas within this 
~ 100 kpc “cool core” (CC) should result in a long-lived cascade 
of multiphase clouds collapsing into the galaxy at its center, fu¬ 
eling extreme star formation rates (10 2 - 10 3 M 0 yr -1 ) amid mas¬ 
sive reservoirs (~ 10 12 M 0 ) of cold molecular gas (e.g., review by 
Fabian et al. 1994b). Although brightest cluster galaxies (BCGs) 
embedded in CC clusters do preferentially harbour these supposed 
cooling flow mass sinks, the observed star formation rates and 
cold gas masses are often orders of magnitude below predictions, 
and high resolution X-ray spectroscopy of the intracluster medium 
(ICM, e.g., Sarazin 1986) is only consistent with reduced cooling 
at ~ 10% of the expected classical rates (e.g., review by Peterson 
& Fabian 2006). 

The mechanical dissipation of active galactic nucleus (AGN) 
power is now routinely invoked by theorists and observers as a so¬ 
lution to the problem, as the average associated energy budget for 
groups and clusters is large enough to inhibit or replenish cool¬ 
ing flow radiative losses not only at late epochs (e.g., Birzan et al. 
2004, 2008; Rafferty et al. 2006; Dunn & Fabian 2006; Best et al. 
2006, 2007; Mittal et al. 2009; Dong et al. 2010), but perhaps 
over a significant fraction of cosmic time (e.g., Hlavacek-Larrondo 
et al. 2012; Simpson et al. 2013; McDonald et al. 2013b). The 
paradigm is motivated by strong circumstantial evidence, including 
nearly ubiquitous observations of radio-bright AGN outflows driv¬ 
ing shocks and excavating kpc-scale buoyant cavities in the ambient 
X-ray gas, acting as lower-limit calorimeters to the often extreme 
(<^ 10 46 ergs sec -1 ) AGN kinetic energy input (e.g., reviews by Mc¬ 
Namara & Nulsen 2007, 2012; Sun 2012). Yet amid panoramic 
supporting evidence (reviewed by Fabian 2012), the physics that 
govern the spatial distribution and thermal coupling of AGN me¬ 
chanical energy to the multiphase (10 -10 7 K) gaseous environment 
remain poorly understood, and cooling flow alternatives invoking 
(e.g.) wet mergers, thermal conduction, and evaporation have been 
a persistent matter of debate (e.g., Bregman & David 1988; Sparks 
1992, 1997; Sparks et al. 1989, 2009, 2012; Fabian et al. 1994a; 
Soker 2003; Voit et al. 2008; Voit 2011; Smith et al. 2013; Canning 
et al. 2015; Voit & Donahue 2015). 

Although often invoked exclusively as a star formation 
quenching mechanism, observations have long demonstrated that 
AGN mechanical feedback does not completely offset radiative 
losses or establish an impermeable “entropy floor”, instead permit¬ 
ting residual cooling either at constant low (~ 10%) rates, (e.g., 
Tremblay et al. 2012a, b), or in elevated episodes as the AGN varies 
in power (e.g., O’Dea et al. 2010; Tremblay 2011). Relative to field 
galaxies or those in non-cool core clusters, BCGs in cool cores pref¬ 
erentially harbour radio sources and kpc-scale filamentary forbid¬ 
den and Baimer emission line nebulae amid 10 9 - 10 11 M 0 repos¬ 
itories of vibrationally excited and cold molecular gas (Heckman 
1981; Hu et al. 1985; Baum 1987; Heckman et al. 1989; Burns 
1990; Jaffe & Bremer 1997; Donahue et al. 2000; Edge 2001; Edge 
& Frayer 2003; Salome & Combes 2003; McNamara et al. 2004; 
Egami et al. 2006; Salome et al. 2006; Edwards et al. 2007; von der 
Linden et al. 2007; Salome et al. 2011; Wilman et al. 2009; Edge 
et al. 2010a,b; McNamara et al. 2014; Russell et al. 2014). Low to 
moderate levels (^ 1- ;> 10 M 0 yr -1 ) of star formation appear to be 
ongoing amid these mysteriously dusty (Quillen et al. 2008; O’Dea 
et al. 2008; Edge et al. 2010a,b; Mittal et al. 2011; Tremblay et al. 


2012a; Rawle et al. 2012), PAH-rich (Donahue et al. 2011) cold 
reservoirs on 50 kpc scales in clumpy and filamentary distribu¬ 
tions (e.g., Johnstone et al. 1987; McNamara et al. 2004; O’Dea 
et al. 2004, 2008; Rafferty et al. 2006, 2008; O’Dea et al. 2010; 
McDonald et al. 201 lb; Tremblay et al. 2014). The ionisation states 
of the nebulae have been a mystery for three decades, and debate 
continues over the roles played by stellar photoionization, shocks, 
thermal conduction, mixing, and cosmic ray heating (e.g., Voit & 
Donahue 1997; Ferland et al. 2009; Sparks et al. 2009,2012; O’Dea 
et al. 2010; McDonald et al. 2010, 2011a; Tremblay 2011; Fabian 
et al. 2011b; Mittal et al. 2011; Oonk et al. 2011; Johnstone et al. 
2012). Whatever the case, there is strong evidence that young stars 
might play an important (although not exclusive) role in dictating 
the physics of both the warm (~ 10 4 K) and cold (~ 100 K) phases 
of the filaments (Voit & Donahue 1997; O’Dea et al. 2010; McDon¬ 
ald et al. 2010, 2011a; Canning et al. 2010, 2014). 

Recent work on star formation in CC BCGs has demonstrated 
its efficacy as an observable tracer for otherwise unobservable 
physical processes regulating the heating and cooling balance in 
hot atmospheres. While low in general and effectively zero in some 
cases, the observed star formation rates in CC BCGs are sometimes 
high enough (^ 100 M 0 yr -1 ) to match condensation rates from the 
X-ray halo (O’Dea et al. 2008). Emergent work at higher redshift 
has shown that the long-ago-predicted classical cooling flows may 
exist after all, forming stars at many hundreds of solar masses per 
year (i.e. the Phoenix cluster, McDonald et al. 2012, 2013a, 2014; 
see also work on Abell 1068 by McNamara et al. 2004; Wise et al. 
2004). Cooling flows may begin to form stars when the central en¬ 
tropy or cooling time drops below a critical threshold (e.g., Voit & 
Donahue 2005; Rafferty et al. 2008; Cavagnolo et al. 2008; Guo 
& Mathews 2013; Voit et al. 2008, 2015a), or when the ratio of 
cooling-to-dynamical times permits a self-regulating “rain” of ther¬ 
mally unstable, spatially inhomogeneous clouds condensing from 
the hot atmosphere (Sharma et al. 2012b; McCourt et al. 2012; Gas- 
pari et al. 2012; Li & Bryan 2014a,b; Li et al. 2015; Brighenti et al. 
2015; Voit & Donahue 2015; Voit et al. 2015a). There is also some 
observational evidence for enhanced cooling in spatially confined 
“cooling channels” where AGN heating may be locally inefficient 
(see, e.g., evidence for enhanced cooling in regions perpendicular 
to the projected cavity/radio “heating axis” in Perseus and Abell 
2597; Lim et al. 2008; Tremblay et al. 2012a). 

Direct observations of young stars in BCGs can test predic¬ 
tions of these various models. To that end, this paper presents a 
morphological analysis of new and archival Hubble Space Tele¬ 
scope ( HST) far-ultraviolet (FUV) continuum images of young, 
massive stars in 16 low-redshift (z < 0.29) CC BCGs. X-ray, Lya, 
Ha, broadband optical, and radio data are also leveraged to cre¬ 
ate an “atlas” of star formation locales relative to the ambient hot 
(;> 10 7 K) and warm ionised (~ 10 4 K) gas phases, as well as the 
old stellar population and radio-bright AGN outflows. In Section 
2 we discuss the sample selection, observations, and data reduc¬ 
tion. Our results are presented in Section 3, discussed in Section 
4, and summarised in Section 5. An appendix contains additional 
multiwavelength overlay figures for all sources in our sample. We 
will frequently abbreviate target names in an obvious manner (i.e. 
Abell 2597 is written as A2597, etc.). Unless otherwise noted, we 
use the names of the parent clusters to refer to their central BCGs 
(i.e. “Perseus” refers to its brightest cluster galaxy, NGC 1275). 
Cosmology dependent physical quantities quoted in this paper as¬ 
sume a flat ACDM model wherein Hq = 10 km s -1 Mpc -1 , Qm = 0.3, 
and Oa = 0.7. Errors are quoted at the la level, unless otherwise 
noted. 
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Table 1. Basic information on the 16 low-redshift cool core brightest cluster galaxies that make up our sample. (1) Source name (note that while we list the 
most commonly used name for the cluster, the actual target studied in this paper is the central brightest cluster galaxy of the named cluster); (2) non-exhaustive 
list of other commonly used names for the cluster, central brightest cluster galaxy, or central radio source; (3) right ascension and (4) declination for the J2000.0 
epoch; (5) redshift (z) as listed in the NASA/IPAC Extragalactic database (NED); (6) the number of kiloparsecs (kpc) that correspond to one arcsecond at the 
given redshift in our assumed cosmology (Hq = 70 km s -1 Mpc -1 ); (7) the associated figure number(s) where the multiwavelength data for the listed BCG can 
be viewed. The FUV continuum images for the entire sample can be viewed at a glance in Fig. 1. 


Source Name 
(1) 

Associated Name(s) 

(2) 

R.A. 

(J2000) 

(3) 

Dec. 

(J2000) 

(4) 

Redshift 

(z) 

(5) 

kpc/arcsec 

(6) 

Shown in 

Figure # 

(7) 

Abell 11 


OOh 12m 44.8s 

-16° 26' 19" 

0.1660 

2.81 

1,2, Al 

Abell 1068 

ZwCl 1037.6+4013 

lOh 40m 47.1s 

+39° 57' 19" 

0.1375 

2.40 

1,7, A2 

Abell 1664 

RX J1303.7-2414 

13h 03m 41.8s 

-24° 13' 06" 

0.1283 

2.27 

1,2, 8, A3 

Abell 1795 

ZwCl 1346.9+2655 

13h 49m 00.5s 

+26° 35' 07" 

0.0625 

1.19 

1,3, 10, A4 

Abell 1835 

ZwCl 1358.5+0305 

14h 01m 02.0s 

+02° 51' 32" 

0.2532 

3.91 

1,2, 7, 11, A5 

Abell 2199 

NGC 6166, 3C 338 

16h 28m 38.5s 

+39° 33' 06" 

0.0302 

0.60 

1,9, A6 

Abell 2597 

PKS 2322-122 

23h 25m 18.0s 

-12° 06' 30" 

0.0821 

1.53 

1,3, 8,2, 10, 11, A7 

Centaurus 

NGC 4696, Abell 3526, PKS 1245-41 

12h 48m 49.2s 

+41° 18' 39" 

0.0099 

0.20 

1,9, 11, A8 

Hydra A 

Abell 780, 3C218 

09h 18m 05.7s 

-12° 05' 44" 

0.0549 

1.05 

1,3,8, 10, 11, 12, A9 

Perseus 

NGC 1275, Abell 426, 3C 84 

03h 19m 48.1s 

+41° 30'42" 

0.0176 

0.35 

1,3, 11, A10 

PKS 0745-191 


07h 47m 31.3s 

-19° 17' 40" 

0.1028 

1.89 

1, 11, All 

RX J1504.1-0248 


15h 04m 07.5s 

-02° 48' 16" 

0.2153 

3.46 

1,7, A12 

RX J2129.6+0005 


21h 29m 37.9s 

+00° 05' 39" 

0.2350 

3.70 

1,2, A13 

ZwCl 0348 

ZwCl 0104.4+0048 

Olh 06m 58.0s 

+01° 04' 01" 

0.2545 

3.93 

1,2, A14 

ZwCl 3146 

ZwCl 1021.0+0426 

lOh 23m 39.6s 

+04° 11' 10" 

0.2906 

4.32 

1,2, A15 

ZwCl 8193 

ZwCl 1715.5+4229 

17h 17m 19.0s 

+42° 26' 57" 

0.1829 

3.04 

1,2, A16 


2 SAMPLE, OBSERVATIONS, AND DATA REDUCTION 
2.1 Sample selection 

The 16 low redshift (z < 0.3) CC BCGs that make up our sam¬ 
ple are listed in Table 1. All are well studied in the literature, and 
enjoy nearly complete cross-spectrum (radio through X-ray) data 
coverage from many ground- and space-based facilities, including 
the Chandra X-ray Observatory , Hubble Space Telescope , Spitzer 
Space Telescope , and Herschel Space Observatory. Eleven of these 
targets constitute the Herschel cool core clusters Open Time Key 
Project sample of A. Edge and collaborators (Edge et al. 2010a, b; 
Mittal et al. 2011,2012; Tremblay et al. 2012a, b; Rawle et al. 2012; 
Hamer et al. 2014), selected to span a wide range of both cool¬ 
ing flow and BCG physical properties. The remaining five targets 
are from the non-overlapping sample of O’Dea et al. (2010), se¬ 
lected from the Quillen et al. (2008) sample on the basis of elevated 
infrared-estimated star formation rates. 

Although biased, our sample spans decades-wide ranges of X- 
ray mass deposition and star formation rates, Baimer and forbidden 
line luminosities, as well as AGN, radio source, and X-ray cavity 
power (including sources that lack detected cavities). Its constituent 
galaxies therefore occupy unique milestones in the supposed ICM 
cooling and AGN heating feedback cycle over the last ~ 3 Gyr of 
cosmic history (redshifts 0.0099 < z < 0.2906), including sources 
with high and low star formation rates, strong and weak AGN feed¬ 
back signatures, as well as many intermediate locales between these 
extremes. A non-exhaustive summary of these various properties 
can be found in Table 2. 


2.2 New Observations 

Along with more than 50 multiwavelength archival observations, 
this paper presents five new FUV continuum and two new broad¬ 
band optical HST observations, all of which are summarised in Ta¬ 
ble 3. Although much of this data has been re-reduced in a homo¬ 
geneous way for this analysis (see §2.3), we refer the reader to the 


references listed in the rightmost column of Table 3 for more ob¬ 
servational details pertaining to the archival data. 

The new HST FUV and optical images we present were ob¬ 
tained in Cycle 19 as part of General Observer program 12220 
(PI: R. Mittal). The line-free FUV continuum data were obtained 
with the Solar Blind Channel (SBC) MAMA detector of the Ad¬ 
vanced Camera for Surveys (ACS, Clampin et al. 2004). Total ex¬ 
posure times for each target were roughly ~ 2700 sec (roughly one 
HST orbit minus overheads), and the observations were carried out 
with a standard three-point dither pattern. Depending on target red¬ 
shift, we used the F140LP, F150LP, and F165LP long-pass filters 
with pivot wavelengths of 1527 A, 1611 A, and 1758 A, respec¬ 
tively. These filter choices ensured that Lya emission did not fall 
within their bandpasses, which have similar red cutoff wavelengths 
of ~ 2000 A, but different blue cutoff (or minimum) wavelengths 
of 1370 A, 1470 A, and 1650 A, respectively. The plate scale for 
the SBC is 0"034 x 0"030 pixel -1 , and the detector field of view is 
34"6 x 30."8. Although we will not discuss the archival FUV ob¬ 
servations for all other sources, we note that the observing strategy 
employed for those datasets was very similar (if not identical) to 
that used for our new FUV observations. 


To fill a data coverage gap, we also obtained two new line- 
free optical images of Hydra A and RX J1504.1-0248 (hereafter 
R1504) using the UVIS channel of HSTs Wide Field Camera 3 
(WFC3, Dressel 2012). As for the FUV images, a three-point dither 
pattern was used over a ~ 2500 sec total exposure time for each 
target. The F814W and F689M filters with pivot wavelengths of 
8024 A and 6876 A were used for Hydra A and R1504, respec¬ 
tively. The widths of these passbands (1536 A and 683 A) forbid 
optical line contamination (though note that some archival optical 
images we use for other sources do contain optical line emission 
like Ha+[N II]). More details on these new HST observations can 
be found in Mittal et al. (2015, submitted). 
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Table 2. A summary of physical properties of the BCGs in our sample, including their surrounding ~ 100 kpc-scale environment. References for the quantities 
presented here can be found below this table. (1) Target name; (2) Infrared-estimated star formation rate (SFR), typically from Spitzer or Herschel. SFRs 
flagged with a f symbol may suffer from a non-negligible Type II AGN contribution to the IR luminosity, in which case the IR-estimated SFR may be 
somewhat overestimated. Note that SFR estimates vary greatly depending on the method, model, and waveband used. We demonstrate this in column (3), 
which shows an average of all published star formation rates for each source (see Mittal et al. 2015 for the specific values used in these averages); (4) Lowest 
and highest X-ray mass deposition rates that have been published for the listed source. Most of these come from Chandra upper limits, rather than more 
reliable XMM Newton Reflection Grating Spectrometer (RGS) data; (5) Cold molecular hydrogen gas mass in units of x 10 9 Mq \ (6) Central ICM entropy in 
units of keV cm 2 ; (7) 1.4 GHz radio luminosity, based on NVSS or Parkes L-band flux densities; (8) Jet mechanical power roughly estimated from (7) using 
the scaling relation of Cavagnolo et al. (2010); (9) X-ray cavity power estimated from Chandra observations. 


ICM COOLING PROXIES AGN HEATING PROXIES 


Source Name 

(1) 

SFR (IR est.) 
(M q yr -1 ) 
(2) 

Published SFRs 
(M 0 yr 1 ) 

(3) 

Mcool 
(Mq yr -1 ) 

(4) 

Mh 2 

(XlO 9 Mq) 

(5) 

Ko 

(keV cm 2 ) 
(6) 

Pi .4 GHz 

(x 10 40 ergs sec -1 ) 
(7) 

~ Pjet (CIO scaling) 
(x 10 44 ergs sec -1 ) 
(8) 

Pcavity (X ray) 

(x 10 44 ergs sec -1 ) 
(9) 

Abell 11 

35 (a) 

35 


1.1 


9.57 ±0.23 

4.42 


Abell 1068 

188 (a) t 

80 ±60 

40-150 

43 

72 

1.34±0.03 

1.01 


Abell 1664 

15 (a) 

16±7 


23 

14.4 

2.13 ±0.06 

1.43 

0.7 ±0.3 

Abell 1795 

8® 

8 ± 8 

1-21 

6 

19 

11.73 ±0.32 

5.15 

!6.0! 2 ; 3 

Abell 1835 

138 <c) 

119 ± 83 


90 

11.4 

7.85 ±0.16 

3.81 

18.o! 19 0 ° 

Abell 2199 

0.6® 

0.2±0.1 

0-3 

1.4 

13.3 

10.45 ±0.33 

4.72 

2 7 +2 - 5 

Z -0.6 

Abell 2597 

5® 

5 ± 8 

20-40 

1.8 

10.6 

42.04 ± 1.11 

13.42 

1 9+lo 

1 ‘ -0.8 

Centaurus 

0.2® 

0.2±0.1 

2.6-2.9 

< 1 

2.25 

1.16±0.03* 

0.91 

0.08+fj* 

Hydra A 

4(b) 

8 ±7 

11-21 

3.2 

13.3 

395.46± 11.37 

72.08 

6.5 ±0.5 

Perseus 

24® 

30 ±23 

12-29 

8.5 

19.4 

21.54 ±0.63 

8.13 


PKS 0745-191 

17® 

70 ±94 

80-260 

4.5 

12.4 

85.53 ±2.58 

22.86 

17.0®;® 

RX J1504.1-0248 

140® 

237 ±92 


10 

13.1 

11.05 ±0.34 

4.93 


RX J2129.6+0005 

13® 

9±5 



21 

5.70±0.19 

3.00 


ZwCl 0348 

52® 

52 




0.44 ±0.01 

0.44 


ZwCl 3146 

67® 

67 ±59 

420-780 

80 

11.4 

2.60±0.12 

1.66 

58.0^;0 

ZwCl 8193 

59® 

59 




16.75 ±0.38 

6.73 



References. — (1) IR-estimated star formation rates are adopted from: (a) O’Dea et al. (2008); (b) Hoffer et al. (2012); (c) McNamara et al. (2006); (d) 
Donahue et al. (2007); (e) Mittal et al. (2012); and (f) Ogrean et al. (2010). (3) Mittal et al. (2015); (4 & 9) X-ray Mass deposition rates and cavity powers are 
collected from Birzan et al. (2004); Dunn & Fabian (2006); Rafferty et al. (2006); Wise et al. (2007); Tremblay et al. (2012a); Kirkpatrick et al. (2009). (5) Cold 
molecular gas masses are adopted from Edge (2001); Edge & Frayer (2003); Salome & Combes (2003); Tremblay et al. (2012a). (6) Central ICM entropy Kq 
is adopted from the main table of the ACCEPT sample, see e.g. Cavagnolo et al. (2009); (7) X-corrected 1.4 GHz luminosities are based on flux densities from 
NVSS (Condon et al. 1998), except in the case of Centaurus, which uses the 1.4 GHz flux density from the Parkes Radio Telescope. The flux-to-luminosity 
conversion is given in Section 2.3; (8) Cavagnolo et al. (2010). 


2.3 Data Reduction 

All new or archival HST FUV and optical data used in this anal¬ 
ysis were retrieved from either the Mikulski Archive for Space 
Telescopes (MAST 1 ) or the Hubble Legacy Archive (HLA 2 ), 
and MAST products were reduced using the standard on-the- 
fly recalibration (OTFR) pipelines. Chandra X-ray observa¬ 
tions were obtained as level one products from the Chandra 
Data Archive 3 . Exposures were reduced, reprojected, exposure- 
corrected, and merged using using the standard CIAO (Fruscione 
et al. 2006) v4.5 scripts (chandra_repro, repro ject_obs, 
flux_obs) with v4.5.5.1 of the calibration database. Finally, 
while many high resolution Very Large Array (VLA) radio maps 
were kindly provided by colleagues, some raw datasets (for A1068, 
RX J1504, and PKS 0745) had to be obtained from the NRAO 
Archive 4 . The NRAO AIPS 5 package was used for (self-) calibra¬ 
tion, imaging, and deconvolution of these data. 

This paper also presents some datasets that have not been 


1 http://archive.stsci.edu/ 

2 http://hla.stsci.edu/ 

3 http://asc.harvard.edu/cda/ 

4 https://archive.nrao.edu/ 

5 http://www.aips.nrao.edu/ 


re-reduced for this analysis. Data reduction and continuum- 
subtraction details for the HST ACS/SBC Ly a images shown in 
Fig. 2 can be found in O’Dea et al. (2010). Reduction of the 
Maryland-Magellan Tunable Filter (MMTF) narrowband Ha maps 
shown in Fig. 3 is described in McDonald et al. (2010). The 1.4 
Msec Chandra X-ray map of Perseus (shown in Fig. A10) is dis¬ 
cussed at length in Fabian et al. (2011a). L-band radio luminosi¬ 
ties quoted in Table 2 use flux densities from the NRAO VLA 
Sky Survey (NVSS, Condon et al. 1998) assuming the relation 
Pi a GHz = 47rL>£S' I y 0 zzo(l +z) a_1 , where Dl is the Luminosity Dis¬ 
tance to the source, S U{) is the 1.4 GHz radio flux density integrated 
over the source area, vo = 1.4 GHz is the frequency of the obser¬ 
vation, and a is the radio spectral index used in the (negligible) 
^-correction, assumed here to be a = 0.8 if S u oc u~ a (these lu¬ 
minosities are insensitive to choice of a given the narrow and low 
redshift range of our targets). 

All images were spatially aligned using IRAF shifting and reg¬ 
istration tasks. To aid viewing of certain X-ray or optical morpho¬ 
logical features, in many cases we show unsharp masks wherein the 
“smooth” X-ray or optical light has been subtracted from the sur¬ 
face brightness map, highlighting residual edge structures. X-ray 
unsharp masks were made in the CIAO environment by gaussian 
smoothing exposure-corrected maps with both small and large ker¬ 
nel sizes. The heavily smoothed map was then subtracted from the 
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Table 3. A summary of the new and archival observations used in this analysis. Those targets for which new FUV continuum or optical data are presented are 
highlighted in boldface. Where applicable, a reference is given to the earliest publication in which the data were first directly analysed. 


Source Name 

A Regime 

Facility 

Inst. / Mode 

Exp. Time / RMS noise 

Obs. / Prog. ID 

Reference / Comment 

Abell 11 

FUV Cont. 

HST 

ACS/SBC F150LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F125LP 

1170 sec 

11230 

O’Deaetal. (2010) 


Optical 

HST 

WFPC2 F606W 

800 sec 

8719 

O’Dea et al. (2010) 


X-ray 

ROSAT 




All Sky Survey 


8.46+1.46 GHz radio 

VLA 

A,B array 

74/iiJy 

AB0878, AL0578 

O’Dea etal. (2010) 

Abell 1068 

FUV Cont. 

HST 

ACS/SBC F150LP 

2766 sec 

12220 

New; Mittal et al. (2015) 


Optical 

HST 

WFPC2 F606W 

600 sec 

8301 



X-ray 

Chandra 

ACIS-S 

26.8 ksec 

1652 

McNamara et al. (2004); Wise et al. (2004) 


8.46 GHz radio 

VLA 

A array 


AE0117 


Abell 1664 

FUV Cont. 

HST 

ACS/SBC F150LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F125LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Optical 

HST 

WFPC2 F606W 

1800 sec 

11230 

O’Dea et al. (2010) 


X-ray 

Chandra 

ACIS-S 

36.6 ksec 

7901 

Kirkpatrick et al. (2009) 


4.86 GHz radio 

VLA 

C array 

100/LzJy 

AE0099 

O’Dea etal. (2010) 

Abell 1795 

FUV Cont. 

HST 

ACS/SBC F140LP 

2394 sec 

11980 



Optical 

HST 

WFPC2 F555W 

1600 sec 

5212 



Ha Narrow Band 

Baade 6.5m 

IMACS / MMTF 

1200 sec 


McDonald & Veilleux (2009) 


X-ray 

Chandra 

ACIS-S 

30 ksec 

10900 etc. 



8.4 GHz radio 

VLA 

A,C,A/D arrays 

18 hrs 

AG0273 

Ge & Owen (1993) 

Abell 1835 

FUV Cont. 

HST 

ACS/SBC F165LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F140LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Optical 

HST 

WFPC2 F702W 

7500 sec 

8249 



X-ray 

Chandra 

ACIS-I 

117.9 ksec 

6880 






36.3 ksec 

6881 






39.5 ksec 

7370 



4.76 GHz radio 

VLA 

A,C arrays 

2 hrs / 47 p Jy 

AT0211 

Govoni et al. (2009) 

Abell 2199 

FUV Cont. 

HST 

ACS/SBC F140LP 

2161 sec 

12220 

New; Mittal et al. (2015) 


Optical 

HST 

WFPC2 F555W 

5200 sec 

7265 



X-ray 

Chandra 

ACIS-I 

120 ksec 

10748 etc. 

Nulsen et al. (2013) 


5 GHz radio 

VLA 

B,C,D arrays 

6 hrs 

AG0269 

Ge & Owen (1994) 

Abell 2597 

FUV Cont. 

HST 

ACS/SBC F150LP 

8141 sec 

11131 

Oonk et al. (2010); Tremblay et al. (2012a) 


Lya 

HST 

STIS F25SRF2 

1000 sec 

8107 

O’Dea et al. (2004); Tremblay et al. (2012a) 


Optical & Ha 

HST 

WFPC2 F702W 

2100 sec 

6228 

Holtzman et al. (1996) 


Ha Narrow Band 

Baade 6.5m 

IMACS / MMTF 

1200 sec 


McDonald et al. (201 la,b) 


X-ray 

Chandra 

ACIS-S 

39.8 ksec 

922 

McNamara et al. (2001); Tremblay et al. (2012a,b) 





112 ksec 

6934, 7329 

Tremblay et al. (2012a,b) 


8.4 GHz radio 

VLA 

A array 

15 min 

AR279 

Sarazin et al. (1995) 


330 MHz radio 

VLA 

A array 

3 hrs 

AC 674 

Clarke et al. (2005) 

Centaurus 

FUV Cont. 

HST 

ACS/SBC F150LP 

1780 sec 

11681 

Mittal et al. (2011) 


Optical 

HST 

ACSAVFC F814W 

8060 sec 

9427 

Harris et al. (2006) 




ACSAVFC F435W 

8654 sec 

9427 

Harris et al. (2006) 




WFC3 F160W 

392 sec 

11219 

Baldi et al. (2010) 


X-ray 

Chandra 

ACIS-S 

200 ksec 

504,5310,4954,4955 

Fabian et al. (2005) 


1.46 GHz radio 

VLA 

A,B/A 

1.5 hr 

AT211 

Taylor et al. (2002) 

Hydra A 

FUV Cont. 

HST 

ACS/SBC F140LP 

2709 sec 

12220 

New; Mittal et al. (2015) 


Optical 

HST 

ACSAVFC F814W 

2367 sec 

12220 

New; Mittal et al. (2015) 


Ha Narrow Band 

Baade 6.5m 

IMACS / MMTF 

1200 sec 


McDonald et al. (2010, 201 lb) 


X-ray 

Chandra 

ACIS-S 

196 ksec 

4969,4970 

Nulsen et al. (2005) 


4.6 GHz radio 

VLA 

A,A/B,B,C,D 

9 hrs 

AL0032 

Taylor et al. (1990) 

Perseus 

FUV Cont. 

HST 

ACS/SBC F140LP 

2552 sec 

11207 

Fabian et al. (2008) 




STIS F25SRF2 

1000 sec 

8107 

Baum et al. (2005) 


Optical & Ha 

HST 

ACSAVFC F625W 

4962 sec 

10546 

Fabian et al. (2008) 


Ha Narrow Band 

KPNO WIYN 3.5m 

S2kB CCD / KP1495 

3200 sec 


Conselice et al. (2001) 


X-ray 

Chandra 

ACIS-S & -I 

1.4 Msec 

11713 etc. 

Fabian et al. (2011a) 


1.4 GHz radio 

VLA 

C,D,C/D arrays 

3 min 

AT149A 

Condon et al. (1996) 

PKS 0745-191 

FUV Cont. 

HST 

ACS/SBC F140LP 

2715 sec 

12220 

New; Mittal et al. (2015) 


Optical 

HST 

WFPC2 F814W 

1200 sec 

7337 

Sand et al. (2005) 


X-ray 

Chandra 

ACIS-S 

50 ksec 

508,2427 

Fabian (1999) 


8.4 GHz radio 

VLA 

A array 

2 hr 

BT024 

Taylor et al. (1994) 

RX J1504.1-0248 

FUV Cont. 

HST 

ACS/SBC F165LP 

2700 sec 

12220 

New; Mittal et al. (2015) 


Optical 

HST 

WFC3/UVIS F689M 

2637 sec 

12220 

New; Mittal et al. (2015) 


X-ray 

Chandra 

ACIS-I 

40 ksec 

5793 

Bohringer et al. (2005); Ogrean et al. (2010) 


8.46 GHz radio 

VLA 

A array 

3 hr 

AB1161 


RX J2129.6+0005 

FUV Cont. 

HST 

ACS/SBC F165LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F140LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Optical 

HST 

WFPC2 F606W 

1000 sec 

8301 

Donahue et al. (2007) 


X-ray 

Chandra 

ACIS-I 

30 ksec 

9370 



8.46 GHz radio 

VLA 

A array 

50/nJy 

AE117 

O’Dea etal. (2010) 

ZwCl 0348 

FUV Cont. 

HST 

ACS/SBC F165LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F140LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Optical 

HST 

WFPC2 F606W 

700 sec 

11312 

Smith et al. (2010) 


X-ray 

Chandra 

ACIS-S 

50 ksec 

10465 



4.86 GHz radio 

VLA 

A/B array 

66/nJy 

AK359 

O’Dea etal. (2010) 

ZwCl 3146 

FUV Cont. 

HST 

ACS/SBC F165LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F140LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Optical 

HST 

WFPC2 F606W 

1000 sec 

8301 

Donahue et al. (2007) 


X-ray 

Chandra 

ACIS-I 

90 ksec 

909,9371 

Boschin (2002) 


4.86 GHz radio 

VLA 

A/B array 

50/nJy 

ACTST 

O’Dea etal. (2010) 

ZwCl 8193 

FUV Cont. 

HST 

ACS/SBC F150LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Lya 

HST 

ACS/SBC F140LP 

1170 sec 

11230 

O’Dea etal. (2010) 


Optical 

HST 

WFPC2 F606W 

1900 sec 

11230 

O’Dea et al. (2010) 


X-ray 

ROSAT 




All Sky Survey 


8.46 GHz radio 

VLA 

A array 

180/Uy 

AE117 

O’Dea etal. (2010) 
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Figure 1 . HST /ACS SBC images of FUV continuum emission associated with young, massive stars in the central ~ 50 kpc of the 16 low-redshift CC BCGs 
in our sample. It is immediately obvious that the star formation in these systems does not occur in monolithic slabs of gas, but rather in highly complex 
filamentary and clumpy distributions. These images are discussed in general in §3. The five sources for which we are showing newly obtained FUV data are 
marked with an asterisk. All panels are rotated such that east is left and north is up. This paper presents multiwavelength data for each of these sources, which 
can be seen by referring to the Figures listed in column (7) of Table 1 . 


lightly smoothed map, and the residual image was normalized by 
the sum of both smoothed maps. Unsharp masks of the FUV and 
optical HST data were made using essentially the same technique 
in the IRAF environment. 


2.4 The SBC red leak and other contaminants 

The ACS Solar Blind Channel suffers from a poorly characterized 
and highly variable red leak (e.g., Ubeda et al. 2012), wherein the 
FUV long pass filters can permit a substantial amount of “red” (i.e. 
optical) interloper flux through the bandpass, contaminating what 
should otherwise be a pure FUV image. The effect is extremely dif¬ 
ficult to correct for in the absence of multi-band UV imaging, as it 
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depends on both time and detector temperature, varying by as much 
as 30% across five consecutive orbits. There is some evidence that 
the effect has decreased since 2008, and current estimates suggest 
that, at worst, it may artificially boost the FUV count rate by about 
10-20% (STScI ACS team, private communication). 

While one must therefore be wary when interpreting SBC 
FUV images of otherwise red, luminous elliptical galaxies, solace 
is found in the fact that the contributors of red leak photons are 
almost exclusively solar- and later type stars. This means that the 
effect can only significantly contaminate by means of a smooth, dif¬ 
fuse, and very faint background whose surface brightness tracks the 
underlying optical isophotes of the host galaxy’s old stellar compo¬ 
nent. We are therefore able to circumvent the issue in this paper 
by quantitatively and qualitatively interpreting only high surface 
brightness, spatially anisotropic FUV-bright clumps and filaments, 
to which the red leak cannot significantly contribute beyond a slight 
increase in count rate that is effectively uniform across such struc¬ 
tures. The same argument applies for highly variable contamination 
from the old stellar “UV upturn” population (see e.g., O’Connell 
1999, for a review). One can therefore be confident that the clumpy 
and filamentary kpc-scale emission ubiquitously seen in our images 
(see Fig. 1) is almost entirely due to young (^10 Myr), massive 
(;> 5 Mq) stars. We nevertheless caution against over-interpretation 
of the smooth, diffuse emission seen in a few of our images (Abell 
2199, for example 6 ), and stress that no quantitative plots or scien¬ 
tific conclusions presented in this paper are based in any way on 
this diffuse emission. 


3 RESULTS 

The FUV continuum images for our full sample are presented in 
Fig. 1. The scales over which FUV emission is detected varies from 
500 pc (Centaurus) to 67 kpc (A1795). Mean and median largest 
angular sizes (L.A.S) are 30 kpc and 33 kpc, respectively. Esti¬ 
mated star formation rates range from effectively zero or <C 0.1 
M q yr” 1 (Centaurus, RX J2129, A2199) to - 150 M 0 yr“‘ (A1068, 
A1835, RX J1504). These and other properties such as cold molec¬ 
ular gas mass, X-ray estimated mass deposition rates, radio source 
power, and X-ray cavity power are summarised for all targets in 
Table 2. 

The FUV morphological analyses in the sections below come 
with an important caveat: FUV emission is highly sensitive to ex¬ 
tinction by dust. The FUV emission that we do detect likely stems 
only from the outermost layers of dense, dusty star forming clouds, 
which are themselves obscured by intervening dust along the line 
of sight. As the FUV is particularly sensitive to young stars less 
than ~ 10 Myr old and more massive than ~ 5 M©, our images 
should be considered instantaneous “snapshots” of ongoing or very 

6 Abell 2199 and Centaurus are the only two sources in our sample for 
which the SBC red leak may make a dominant contribution to observed 
morphology. As a test, we have scaled the V-band optical image of A2199 
by a factor consistent with the known range of reasonable FUV/optical col¬ 
ors (discussed at length in e.g. Oonk et al. 201 1), then subtracted this scaled 
image from the FUV map. We are unable to rule out the possibility that 
all emission seen in the A2199 “FUV” image is a combination of red leak 
from the underlying old stellar component, plus a central FUV point source 
associated with the AGN. The star formation rate in this source may there¬ 
fore be effectively zero. The same is true for Centaurus, however all other 
sources in our sample are dominated by bright FUV clumps and filaments 
whose morphology will be unaffected by the red leak. 


recent unobscured star formation. A detailed treatment of extinc¬ 
tion for a majority subset of our sample is provided by Mittal et al. 
(2015), and will not be discussed here beyond cautioning against 
over-interpretation of observed FUV structures. The clumps and 
filaments we do detect are likely “tips of icebergs”, and smooth, 
diffuse emission may be significantly contaminated by red leak and 
the UV upturn population (as discussed in §2.4). 

It is nevertheless obvious from Fig. 1 that star formation in 
our sample is not occurring amid monolithic slabs of gas. The ob¬ 
served FUV morphologies are instead highly clumpy and filamen¬ 
tary, exhibiting a variety of associations (and sometimes interesting 
77077 -associations) with X-ray, optical, and radio features, as well 
as galaxy properties such as central X-ray entropy and the relative 
strength of AGN feedback signatures. These associations are dis¬ 
cussed in the following sections. 

3.1 Comparison of FUV, Lya, and Ha morphology 

In Fig. 2 we compare a subset of our targets with the continuum- 
subtracted Lya data of O’Dea et al. (2010). Although the mor¬ 
phologies are very similar overall, the Lya is far more extended 
than the underlying FUV continuum. In a simple photon-counting 
exercise using the same FUV data for a subset of our sample, O’Dea 
et al. (2010) demonstrated that the young stellar component traced 
by the FUV continuum can roughly account for the photon budget 
required to photoionise the Lya nebula, although there is unavoid¬ 
ably significant uncertainty in the extinction correction used in this 
analysis. 

That the Lya is far more extended than the FUV continuum 
may be due to a simple sensitivity issue. Lya is far brighter than 
the local FUV continuum, as the average Lya/FUV flux density 
ratio for our sample is roughly ~ 3. There are several examples 
in our sample where outer Lya filaments are detected in a region 
where this Lya/FUV ratio would result in an FUV continuum flux 
that is below the sensitivity limit of the observation (the outer Lya 
filaments in Abell 11 are one example). We are therefore unable 
to rule out the possibility that all Lya emission is cospatial with 
underlying FUV continuum from young stars. Alternatively, it is 
still possible that the Lya is intrinsically more extended than the 
FUV continuum. This would be similarly unsurprising and con¬ 
sistent with many previous observations, as this is attributable to 
Lya’s very high sensitivity to resonant scattering (e.g., Laursen & 
Sommer-Larsen 2007). This can make Lya morphology difficult 
to interpret in a physical sense, though it does serve as an excel¬ 
lent (although “tip-of-the-iceberg”) tracer for neutral hydrogen. We 
leverage these Lya data for this purpose in many of the multiwave¬ 
length comparison figures listed in column (7) of Table 1 . 

In Fig. 3 we compare a subset of our FUV sample with the 
narrowband Ha maps from McDonald & Veilleux (2009); McDon¬ 
ald et al. (2010, 201 la, b) and Conselice et al. (2001). We have 
smoothed the FUV maps with a gaussian in order to degrade their 
spatial resolution to (roughly) match that of the Ha images. Al¬ 
though the FUV and Ha morphologies closely match one another 
(see also e.g., McDonald et al. 201 la), the match is not nearly one- 
to-one. It has been known for many years that some of the Ha fil¬ 
aments in CC BCGs are devoid of a detectable FUV counterpart, 
with Perseus 7 being the most obvious example (e.g., Hatch et al. 

7 One must be wary of confusing star forming filaments in Perseus with 
FUV and blue excess emission from the foreground High Velocity System 
(HVS) that is superimposed along the line of sight. This disrupted galaxy 
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Figure 2. A comparison of Lyct and FUV continuum morphologies for a subset of our sample. In blue we show the HST /ACS SBC continuum-subtracted 
Lyct images from O’Dea et al. 2010, and in orange we show the FUV continuum images. While the general morphologies are very similar, the Lya emission 
is far more extended than the underlying FUV continuum. This result is unsurprising given the sensitivity of Lya emission to resonant scattering. O’Dea 
et al. (2010) demonstrated that the underlying FUV continuum strength was sufficient to fully account for production of the observed Lya emission via stellar 
photoionization. All image pairs are aligned and shown on a common spatial scale, with east left and north up. 


2006; Canning et al. 2010). Moreover, some blue star forming fila¬ 
ments apparently lack cospatial Ha emission (e.g., the “blue loop” 
in Perseus, Fabian et al. 2008; Canning et al. 2010, 2014). Ha 
traces the contemporary star formation rate via the instantaneous 
flux of ionizing photons from the most massive (M* ;> 15M©) O 

is ~ 100 kpc closer in projection and is unrelated to the BCG; see e.g., 
Sanders & Fabian 2007. 


and early B-type stars, while the more heavily extincted UV ex¬ 
cess associated with the photospheres of less massive (M* ;> 5M 0 ) 
young stars can shine long after those most massive stars powering 
the Ha flux are gone. Ha and the FUV therefore sample smaller 
and larger temporal slices (~ 10 6-7 yr vs. ~ 10 8 yr) of the star for¬ 
mation history, respectively. More importantly, many authors have 
demonstrated that the Ha nebulae cannot be heated by star forma¬ 
tion alone (see §1, although this issue is not the focus of our paper). 
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Narrowband Ha Smoothed FUV Continuum 



Figure 3. A comparison of Ha and FUV continuum morphologies for a subset of our sample. In blue we show the narrowband Ho: images from McDonald 
et al. 2010, 201 la,b and Conselice et al. 2001, and in orange we show the FUV continuum images after gaussian smoothing to (approximately) match the 
spatial resolution of the Hen maps. All image pairs are aligned and shown on a common spatial scale, with east left and north up. The morphologies are very 
similar overall, though it is important to note that some Hen filaments lack detected cospatial FUV continuum, and some FUV filaments lack detected cospatial 
Hen emission (though non-detection of course does not necessarily imply absence). Many works have demonstrated that continuum emission from young stars 
can account for a dominant fraction of the ionizing photons needed to power the Ho; nebula (e.g., O’Dea et al. 2010; McDonald et al. 2011a). However, an 
additional ionization mechanism (e.g., shocks, cosmic rays, thermal conduction, etc.) is needed in all cases (see §1). 
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Figure 4. A comparison of projected position angles (PA.) for the FUV, X-ray, optical, IR, and radio isophotal major axes at matching spatial scales. P.A. has 
been measured N through E. If the P.A. of the X-ray, optical, or IR major axis was found to vary strongly as a function of radius, we measured the P.A. of the 
major axis that matched the largest angular size of the FUV emission. Strong alignment between FUV, X-ray, optical, and IR counterparts is observed. Those 
sources that are outliers to any particular trend are labeled in the respective plot. The radio vs. FUV major axis comparison {bottom right ) shows alignment 
only for those sources that exhibit evidence for either jet-triggered star formation or strong dynamical interaction between the radio source and star forming 
gas. The dashed line on the bottom right plot is the one-to-one line. Spearman-rank and Pearson correlation coefficients for these plots are shown in Table 4. 


The slight morphological mismatch between Ha and FUV is there¬ 
fore not necessarily surprising. Fig. 3 only demonstrates that the 
FUV and Ha filaments are roughly cospatial in projection, with 
clear important exceptions. While some authors have shown that 
stars can indeed play a very important role in the ionisation states 
of the filaments (e.g., O’Dea et al. 2010; McDonald et al. 2011b), 
we reiterate that another heat source (acting either alone or in con¬ 
cert with the young stars) is needed (e.g. Voit & Donahue 1997; 
Ferland et al. 2009; Oonk et al. 2011; Fabian et al. 2011a; Mittal 
et al. 2012; Johnstone et al. 2012; Canning et al. 2015) 

3.2 Galaxy-scale position angle alignment of FUV with 
X-ray, optical, and IR major axes 

In Fig. 4 we plot the position angle (P.A., measured N through E) 
of the projected X-ray, optical, IR (3.6 gm), and radio major axis 


Table 4. : The Spearman-rank and Pearson correlation coefficients between 
the FUV postion angle ( X ) and the X-ray, optical, infrared and radio posi¬ 
tion angles (F). These quantities are plotted against one another in Fig. 4. 


Quantity (F) 

Spearman-Rank 

Pearson 

Best Fitting Correlation 

X-ray 

0.81 

0.84 

Y ocX 1±a2 

Optical 

0.73 

0.74 

FocX°- 9 ± 01 

Infarred 

0.72 

0.76 

Y ocX 1±0 ' 2 

Radio 

0 

0.10 



vs. the projected FUV major axis. The FUV major axis was taken 
to be the position angle of the isophote at roughly twice the FUV 
half-light radius in lightly smoothed maps. We then measured the 
X-ray, optical, and IR major axis within the isophote at roughly the 
same radius. Sources that are point like or circular in any of these 
bands have been excluded from that particular plot. 
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Figure 5. The FUV asymmetry index A, as defined by Eq. 1, vs. the FUV 
“clumpiness index” S, defined in Eq. 2. Higher values of A means that the 
FUV surface brightness distribution is more azimuthally asymmetric about 
the galaxy centre. Higher values of S means that the FUV emission features 
more high spatial frequency clumps, and less smooth emission. The two 
indices strongly correlate, such that a galaxy with a high A + S value can 
be considered more spatially anisotropic (in terms of FUV emission) than a 
galaxy with a low A + 5 value. These indices are discussed in §3.4. 


As is evident in Fig. 4, we observe weak-to-strong projected 
alignment between FUV, X-ray, optical, and IR counterparts. Those 
sources that are outliers to any of these trends have been labeled in 
the respective plot. It is possible that the alignment is a reflection 
of the old and young stellar components sharing a common origin 
in the ambient hot gas. Another possibility is that the alignment is 
merely due to the fact that the various components all reside within 
the same gravitational potential, and have had sufficient time to 
dynamically relax, torque toward a common axis, etc. We caution 
against over-interpretation of these apparent projected alignments: 
these are chaotic, messy systems with morphologies that probably 
vary strongly with time. 


3.3 Kiloparsec-scale offsets between FUV and X-ray surface 
brightness peaks 

While FUV and X-ray surface brightness peaks are spatially coin¬ 
cident for the majority of our sample, A1664, A1835, Centaurus, 
PKS 0745, R2129, and Zw3146 show projected offsets of 9 kpc, 
14.8 kpc, 1.3 kpc, 4.65 kpc, 7 kpc, and 11 kpc, respectively. The 
mean and median offsets are are both roughly 8 kpc. Offsets be¬ 
tween the X-ray emission and the optical/IR BCG peak are effec¬ 
tively the same for these targets (as the FUV and optical peaks are 
almost always cospatial, at least within our sample). The offset in 
A1664 has been previously noted by Kirkpatrick et al. (2009) in 
their detailed study of the source. The A1835, R2129, and Zw3146 
offsets were noted by O’Dea et al. (2010). 

The X-ray surface brightness maps for all objects with X- 
ray/FUV photocentroid offsets show large scale (;> 50 kpc) asym¬ 
metries, suggestive of complex gas dynamics in the hot phase. Pho¬ 


tocentroid offsets between the BCG and the cluster X-ray emission 
is a proxy for how close (or how far) the system is to a state of 
dynamical equilibrium, such that the offsets should decrease as the 
cluster evolves (Katayama et al. 2003). Sloshing motions in the X- 
ray gas can nevertheless remain long-lived even after the supposed 
virialisation of the cluster (e.g. Markevitch & Vikhlinin 2007). 
Large optically selected samples of both cool core and non-cool 
core BCGs frequently show median X-ray/BCG offsets of ~ 15 kpc 
(Bildfell et al. 2008; Sanderson et al. 2009; Loubser et al. 2009). 
CC BCGs systematically lie below this median at ~ 10 kpc, which 
is close to the median for those objects (with observed X-ray off¬ 
sets) within our sample (~ 8 kpc). If we include those galaxies in 
our sample that do not show any measurable projected offset (ten 
out of sixteen sources), the sample-wide median and mean pro¬ 
jected offsets are ~ 0 kpc and ~ 3 kpc, respectively. The sources 
exhibiting kpc-scale offsets do not appear to prefer any particular 
galaxy property - instead they inhabit the full range of FUV mor¬ 
phology, star formation rate, radio power, etc. that is spanned by 
our whole sample. 


3.4 Quantifying morphology by Asymmetry and Clumpiness 
indices 

Fig. 1 shows that our sample spans a diverse range of FUV mor¬ 
phologies, including sources that can be described as amorphous 
(e.g., Zw3146), clumpy (e.g., All), point-like (e.g., Centaurus), 
disk-like (Hydra A), and filamentary (e.g., A2597, R1504). There 
is however significant overlap between these classes. For example, 
A2597 could be arbitrarily described as a hybrid of filamentary, 
amorphous, and clumpy structures, illustrating the need for a more 
objective measure of morphology. 

We therefore quantify all projected FUV morphology by 
means of scale- invariant structural indices, as is done frequently 
for galaxies in the literature. We adopt the commonly used 
Concentration-Asymmetry-Smoothness (CAS) system described 
by Conselice (2003), which posits that galaxy morphology can 
be entirely quantified by measuring the concentration of light (C) 
around a photocentric point, the azimuthal asymmetry of light 
about this point (A), and the high spatial frequency smoothness or 
clumpiness ( S ) of that light. The CAS indices are useful in that they 
(a) are independent of any assumption about galaxy light distribu¬ 
tion and (b) correlate with galaxy processes such as star formation, 
mergers, colors, emission line widths, etc. Galaxies of different 
Hubble type appropriately stratify within the optical CAS volume, 
which has been expanded to include other wavelength regimes over 
the years (including for extragalactic FUV imaging, see e.g., Hol- 
werda et al. 2012). 

In our case, there is a risk that any use of a concentration-of- 
light parameter C (typically defined by the ratio of curve-of-growth 
radii containing 80% and 20% of all light) may be significantly 
contaminated by the SBC red leak, for reasons discussed in §2.4. 
We therefore make use of only the asymmetry and clumpiness pa¬ 
rameters A and S , which (even without C) are useful in quantifying 
spatial anisotropy in FUV surface brightness. Following Conselice 
(2003), we compute the asymmetry index A by rotating each FUV 
image by 180° about a central point (discussed below), subtract¬ 
ing this from the original unrotated image, and then summing the 
absolute value intensities from the resulting residual map. The re¬ 
sultant value is then normalized by two times the original galaxy 
flux. Expressing the above more quantitatively, the asymmetry in¬ 
dex is given by 
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where Ie=o and Io=m are the intensity distributions in the original 
and 180° rotated images, respectively. On both the rotated and orig¬ 
inal images, the sum is taken over all pixels within matching re¬ 
gions that encompass all galaxy FUV flux (e.g., an “all the flux you 
see” circular aperture). The resulting value of A depends strongly 
on the pixel about which the image is rotated, as is discussed by 
Conselice (2003). For sample-wide consistency, we have chosen to 
rotate each FUV image around the pixel that is cospatial with both 
the radio core and optical photocentroid of the host galaxy, such 
that our computed A values are at least somewhat related to the 
projected reflection asymmetry of young stars around the AGN. In 
a few cases the radio core was not cospatial with the host galaxy 
optical photocentroid, but this was typically because of a central 
dust lane (Hydra A is one example). Comparison of both A and S 
values with those in other papers is beyond the scope of this work, 
so in the interests of simplicity we use k as a scalar normalization 
to set the range spanned by the A distribution to 0 < A < 1. Sources 
with lower values of A are more azimuthally symmetric about the 
centre of the galaxy than are sources with higher values of A. 

The “dumpiness” parameter S was calculated by summing 
pixel intensities in an unsharp mask of the FUV image, made by 
gaussian smoothing the map with both small and large kernel sizes, 
subtracting the heavily smoothed map from the lightly smoothed 
one, and then normalizing the residuals by the flux in the original 
image. More specifically, S is given by 


S = jx 


(I-Ia)-B 

J 1 


( 2 ) 


where I and I a are the intensity distributions in the lightly and 
heavily smoothed images, respectively. All of these sums are taken 
within a matching area, and all galaxy apertures are made with a 
central hole that intentionally excludes the galaxy nucleus (where 
an FUV point source might artificially weight the S value). B is the 
background intensity distribution within an off-source “sky” aper¬ 
ture. We again normalize by j to set the range of possible S values 
equal to that of A, i.e. 0 < S < 1. A galaxy with with 1 will fea¬ 
ture many high spatial frequency clumps, whereas an object with 
S ~ 0 is very smooth. 

We plot A vs. S in Fig. 5. The two indices strongly correlate, 
such that sources that are more asymmetric also tend to be more 
clumpy. The trend is strong enough that for the remainder of this 
paper we will discuss A and S together by means of a single “FUV 
anisotropy index” A + S, which can inhabit the range from 0 < A + 
S < 2. Galaxies with higher A + S are more clumpy, filamentary, 
and asymmetric, whereas galaxies with lower A +S values are more 
symmetrically amorphous or point-like. 


3.5 Comparison of FUV morphology with redshift, star 
formation rate, and ICM central entropy 

FUV morphology does not exhibit any obvious redshift depen¬ 
dence, despite the very strong redshift-luminosity bias present in 
our sample (which is assembled from flux-limited and therefore 
Malmquist-biased catalogs). We demonstrate this bias in the top- 
left panel of Fig. 6, where we plot the IR-estimated star formation 
rate vs. redshift. Despite the expected strong upward trend in SFR 
with redshift associated with the Malmquist bias, the top right panel 
of Fig. 2 — showing FUV anisotropy index (A +S) vs. redshift — 
is effectively a scatter plot (note that, for each of these plots, we 
“gray out” red-leak contaminated A2199 and Centaurus, so they 


do not give the illusion of correlation). Galaxies at higher redshift 
marginally tend to have a higher FUV anisotropy value, albeit with 
very large scatter. The error bar on this plot reflects the rather large 
range that A+ 5 can inhabit given slightly different choices of pixel 
about which the image is rotated (in the case of A) or smoothing 
lengths used to make the unsharp mask (in the case of S). 

We conclude that there is no evidence for any correlation 
between redshift and morphology in our sample. This is perhaps 
somewhat surprising, because one naturally expects a trend be¬ 
tween redshift and morphology as (1 +z) 4 surface brightness dim¬ 
ming and angular size scaling should make objects look increas¬ 
ingly smooth and symmetric as they approach the resolution limit 
at higher redshifts. To independently test what effect redshift may 
have on perceived morphology in our sample, we have artificially 
redshifted all of our FUV images to one common redshift equal 
to that for our most distant target (Zw3146, z = 0.2906). An IRAF 
script was used to accomplish this, implementing the technique de¬ 
scribed by Giavalisco et al. (1996) (specifically, see their equations 
2-7). While artificially redshifting our targets had only a small ef¬ 
fect on overall asymmetry A, it is clear from this test that redshift 
can have a strong effect on perceived smoothness S. Our lowest red¬ 
shift targets around z = 0.01, for example, suffer a factor of ~ 20 
degradation in spatial resolution, lowering their A value negligibly 
and their S value moderately (depending on choice of smoothing 
scale lengths). It is therefore possible that our high redshift targets 
are intrinsically far more clumpy than they appear. We therefore 
note that, while the top right panel of Fig. 6 provides no evidence 
for a correlation between A + S and redshift, it cannot be used to 
rule it out. 

We plot IR-estimated star formation rate vs. FUV anisotropy 
index in the lower left panel of Fig. 6, finding no correlation. We 
do observe a weak upward trend of central ICM entropy with FUV 
anisotropy index, as is evident from the lower right panel of Fig. 6. 
Here, entropy S (in units of keV cm 2 ) is defined as S = kTn ^ 2//3 , 
where k is the Boltzmann constant, T is the gas temperature, and n e 
the electron density. Central entropy values have been adopted from 
the ACCEPT sample (Donahue et al. 2006; Cavagnolo et al. 2009). 
We again caution against over-interpretation here, particularly be¬ 
cause calculations of ICM central entropy from the X-ray data can 
be problematic and strongly tied to data quality (Panagoulia et al. 
2014). The plot merely demonstrates that sources with higher cen¬ 
tral entropy may have more spatially anisotropic star formation, at 
least within our sample. 


3.6 Star forming filaments aligned with radio jets and lobes 

Fig. 7 through Fig. 9 shows selected X-ray, FUV, Lyce, and radio 
overlay figures for a subset of our sample. The three figures are 
presented in order of highest to lowest star formation rate, respec¬ 
tively. In these panels one will find several clear spatial correlations 
between FUV emission, radio emission, and/or X-ray emission. We 
discuss these correlations in the next two sections. 

In Fig. 10 we highlight four examples of strong morphological 
alignment between FUV continuum / line emission and radio jets 
or lobes (shown in red contours). These include A1795, Hydra A, 
and A2597. While A1795 and A2597 are known and well-studied 
examples (e.g., O’Dea et al. 2004), such alignment has not previ¬ 
ously been noted for Hydra A. Moreover, Fig. All shows some 
evidence of alignment in PKS 0745, whose “spike-like” FUV fila¬ 
ment is aligned with the axis about which the radio source appears 
to kink or fold over. 
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Figure 6. (top left ) Infrared-estimated star formation rate vs. redshift for all targets in our sample. SFRs estimated by other indicators in other wavelength 
regimes tend to be lower than the IR-based rate, which can be considered a rough upper-limit. (top right) “FUV anisotropy index” (A + 5) based on the CAS 
parameters described by Conselice (2003). We define the anisotropy index in §3.4. A2199 and Centaurus have been “grayed out” so as not to give the illusion 
of correlation where there is likely none. These two sources may be highly contaminated by red leak (so their CAS morphology cannot be trusted), and their 
star formation rates are extremely low (perhaps effectively zero). 


There are now too many examples of clumpy /filamentary star 
formation alinging with radio jets and lobes for this cospatiality to 
credibly be pure coincidence or a projection effect (see e.g., Cen 
A — Crockett et al. 2012; Hamer et al. 2015; Santoro et al. 2015; 
Minkowski’s Object / NGC 541 — van Breugel et al. 1985; 3C 285 
— van Breugel & Dey 1993; 4C 41.47; Bicknell et al. 2000; see 
also 3C 305, 3C 321, 3C 171, and 3C 277.3). These filaments may 
have been dynamically entrained, uplifted, or swept aside by the 
radio source, or may have formed in situ along the working sur¬ 
face of the radio lobe in an example of positive AGN feedback. 
It has long been predicted that star formation may be triggered by 
shock- induced cloud collapse as the propagating radio plasma en¬ 
trains and displaces cold gas phases (see e.g., the shock/jet-induced 
star formation models by Elmegreen & Elmegreen 1978; Voit 1988; 
De Young 1989; McNamara & O’Connell 1993). Jet-induced star 
formation has for many years been considered as a plausible ex¬ 


planation for the high redshift alignment effect (Rees 1989; Daly 
1990). 

Whatever the case, examples like these demonstrate that, at 
least for some time, star formation can survive (and may indeed 
be triggered by) dynamical interaction with a propagating radio 
source. If jet-triggered star formation is indeed a real effect, and is 
responsible for the alignment seen for the four targets in our sam¬ 
ple, we can roughly estimate whether or not such an effect has a sig¬ 
nificant or negligible impact on the global star formation rate in the 
galaxy. If we take A1795, our most dramatic example, and assume 
that all FUV emission associated with the “P”-shaped filament 
cospatial with the radio lobes is directly induced by propagation 
of the jet, then up to 50% of all star formation in the galaxy could 
be jet-triggered. The upper limit percentages for the other sources 
in our sample are much lower, such that even if jet-triggered star 
formation is indeed real, it probably does not play the dominant 
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Figure 7. The three sources in our sample with the highest star formation rates (:> 100 Mq yr -1 )- The left-hand panels show a wide view of X-ray emission 
cospatial with the BCG and its outskirts. White boxes are used to indicate the FOV of the right-hand panels, which show FUV continuum emission. Various 
contour sets are overlaid, and are labeled appropriately in their respective panels. 
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Figure 8. A selection of sources in our sample that show moderate star formation rates ~ 5 < SFR <15 Mq yr _1 . The left-hand panels show a wide view 
of X-ray emission cospatial with the BCG and its outskirts. White boxes are used to indicate the FOV of the right-hand panels, which show FUV continuum 
emission. Various contour sets are overlaid, and are labeled appropriately in their respective panels. 
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Figure 9. The two sources in our sample in which there is effectively no ongoing star formation or FUV-bright young stars (or <0.1 A/ 0 yr -1 )-The left-hand 
panels show a wide view of X-ray emission cospatial with the BCG and its outskirts. White boxes are used to indicate the FOV of the right-hand panels. For 
Centaurus, we show a B/H -band color map made with HST data, showing the famous dust lane associated with the source. The (almost nonexistent) FUV 
emission in the radio core is shown in red contours. For A2199, the rightmost panel shows FUV continuum emission. The diffuse emission seen in A2199 may 
be an artifact of the ACS SBC red leak. 


role in driving star formation in the galaxy. We estimate that the 
effect may play a role at the few percent level at best. Regardless, 
the apparently competing roles of radio mechanical feedback si¬ 
multaneously quenching and triggering star formation can be rec¬ 
onciled with one another. Even if the propagating radio source does 
not (immediately) inhibit or truncate star formation directly, it may 
still work to starve it of gas for future star formation by excavating 
cavities and driving sound waves in the hot gas, preventing it from 
cooling and forming stars. 


3.7 Spatial correlations and anti-correlations of star forming 
filaments with X-ray cavities 

As we demonstrate in Fig. 11, six sources in our sample possess 
one or more kpc-scale narrow filaments that, in projection, extend 
toward, into, or wrap around the edges of kpc-scale X-ray cavi¬ 
ties. These include Perseus, A2597, Hydra A, PKS 0745, Centau¬ 
rus, and A1835. Two additional sources (A1664 and A1068) show 
weaker evidence (due perhaps to the unavailability of deeper X-ray 
imaging) of the same effect. Most of those filaments that extend 
toward and into cavities are FUV bright and forming stars, while 
Perseus and Centaurus show only dusty, Ha bright filaments that 
lack cospatial FUV continuum. 

This may be evidence for buoyant uplift of the filament by the 
cavity as it rises amid the ICM, as has been discussed (typically 
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Figure 10. The strongest examples of FUV/radio morphological correlation and anticorrelation in our sample. This is perhaps evidence for for (a) star forming 
filaments that have been uplifted, entrained, or swept aside by the propagating radio source or (b) jet-triggered star formation. Note also that PKS 0745 
(Fig. All) shows weaker evidence for a similar alignment. 



Figure 11. Six sources for which there is some evidence of kpc-scale filaments extending in projection towards and around X-ray cavities. FUV contours are 
overlaid in green on X-ray unsharp masks. Centaurus has no discernable star formation, but its dust contours do possess filaments that extend toward cavities. 
We show these dust contours in blue. For Hydra A we show Ho: contours as these better show the faint filament that follows the northern cavity. While not 
shown, FUV emission at low surface brightness in Hydra A is fully cospatial with the Hen filaments shown above. 
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Figure 12. The new HST /-band optical image of the Hydra A BCG and its surrounding environment. On the right-hand panel we overlay both the VLA 
4.6 GHz radio contours and the MMTF narrowband Ha contours. The Ha distribution shows evidence for dynamical interaction with both the radio source 
(note the apparently uplifted filaments northward and southward of the nucleus) as well as the small companion galaxy ~ 10" to the southeast. While this 
companion is unlikely to to provide a substantial cold gas mass to the BCG via merger-driven flow, it may be acting to dynamically “stir” the low entropy gas 
already present in the BCG. The companion is cospatial with a bright knot of both Ha and FUV continuum emission. 


in the context of Perseus) by many previous authors (e.g., Fabian 
2003; Hatch et al. 2006; McDonald et al. 2011b; Canning et al. 
2010, 2014). The spatial associations are certainly compelling, and 
filament uplift by cavities may indeed be an important and even 
common effect. It is, however, unlikely to be the only effect driving 
the morphology of all narrow filaments ubiquitously observed in 
CC BCGs, as there are many examples of filaments with no obvi¬ 
ous association with either a cavity (or radio source, for that mat¬ 
ter). Ever-deeper observations of X-ray cool cores do however tend 
to reveal ever more numerous X-ray cavities, so we cannot neces¬ 
sarily rule out the unlikely possibility that all filaments in CC BCGs 
have at some point been uplifted by a cavity. We find this unlikely, 
however, as filament kinematics (at least for the small sample that 
has been fully mapped with an IFU) are generally inconsistent with 
expectations if they are indeed dragged outwards (see, for exam¬ 
ple, the northern filament in Perseus, which shows a smooth veloc¬ 
ity gradient of a few 100 km sec -1 , consistent with a laminar flow; 
Hatch et al. 2006). Alternatively, those filaments that wrap (in pro¬ 
jection) around X-ray cavities may have formed in situ in the cav¬ 
ity’s compressed shell, though it is entirely unknown whether or not 
direct cooling from the X-ray to molecular phase is even possible 
absent dust (e.g., Fabian et al. 1994b). We expand on this below, in 
our discussion of FUY morphology in the contex of ICM cooling 
and AGN heating. 


4 DISCUSSION 

The exquisitely complex and highly diverse range of far ultravi¬ 
olet morphologies prestented in this paper reflect the dynamical 
response of low entropy gas to a highly energetic, chaotic envi¬ 
ronment. Besides the graviational potential of their host galaxies, 
these star forming nebulae reside amid AGN-driven jets, bubbles, 
sound waves, bulk ICM motions, and stellar feedback, and few 
structures observed in our FUV dataset are likely to be long-lived. 
Hydra A’s rotating, star forming disk (Hamer et al. 2014) reveals 
cold gas that is largely in dynamical equilibrium, though it still fea¬ 


tures narrow filaments that have likely been lifted outward by the 
radio jet (Figs. 10 & A9), or dynamically stirred by a small compan¬ 
ion galaxy (Fig. 12). A1795 (Fig. A4) not only features radio lobes 
frosted with young stars, but a ~ 20 kpc southern tail deposited 
perhaps by a cooling wake that lags behind the BCG (McDonald 
& Veilleux 2009). Centaurus (Fig. A8) features no discernable star 
formation whatsoever, but a spectacular winding dust lane whose 
shape mirrors that of the larger scale X-ray spiral. 

In many ways, conclusions drawn from a small collection of 
highly complex individual galaxies are doomed to ambiguity. If 
each source is so chaotic and time-varying, what can we learn in 
a “big picture” context? The answer may be that our data are snap¬ 
shots of a prototypical cool core BCG at different stages of an AGN 
outburst cycle. Throughout the past ~ 3 Gyr of cosmic history that 
our sample’s redshift range spans, perhaps each galaxy has spent 
(or will spend) some time resembling each of the others as the tug- 
of-war between ICM cooling and AGN heating cycles and varies. 
Indeed, AGN can vary in power and switch on and off over a man¬ 
ifold range of timescales, and the associated balance of AGN heat¬ 
ing and ICM cooling can vary still more. 

Recent numerical work by Sharma et al. (2011); McCourt 
et al. (2012), and Gaspari et al. (2012) has shown that thermal 
instabilities in a cooling flow can produce a multiphase and star 
forming ISM when the ratio of the cooling time t co oi to the local 
gravitational free-fall timescale fir is Cool /Cf 10 (see also Pizzo- 
lato & Soker 2005, 2010). This theoretical framework has since 
been expanded in a series of papers by Voit and collaborators, who 
propose a precipitation-regulated AGN feedback model applicable 
not only to BCGs and giant ellipticals (Voit & Donahue 2015; Voit 
et al. 2015a, b), but perhaps galaxies in general (Voit et al. 2015c, in 
press). In the model, cold clouds precipitate out of the ambient hot 
medium via thermal instability wherever Cool 10fif. Cold chaotic 
accretion (Gaspari et al. 2013, 2014) from the now “raining” am¬ 
bient hot atmosphere boosts black hole feeding to ~ 100 times the 
Bondi rate, powering jets that can stimuluate further precipitation 
by dragging low entropy gas to higher altitudes, where the cooling- 
to-dynamical time ratio will lessen. At the same time, jet heating 
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Figure 13. Here we compare the maximal FUV radius with the X-ray entropy (left) and cooling-to-dynamical time ratios (right) measured at that same cluster¬ 
centric radius. Assuming (a) that the star formation in our sources is indeed powered by a cooling flow and (b) there is no non-detected star formation beyond 
the largest measured FUV radius, the maximal radius within which FUV emission is observed serves as a rough observable tracer of the radial threshold for the 
onset of cooling flow powered star formation. Within our sample, the average maximal radius (and its ±lcr interval) within which FUV continuum emission 
is detected is 14.7 =b 10.14 kpc. At this average radius, the average (and drier) cooling time, entropy, and cooling-to-dynamical time ratio is 1.09 ±0.68 
Gyr, 28.27 ± 12.12 keV cm 2 , and 14.78 ± 8.57, respectively. These average values and their ±lcr intervals are marked by the black solid and dashed lines 
(respectively) on both panels. Blue and red points are used to differentiate between those sources with unresolved radio sources (and higher SFRs) and resolved 
radio sources (and lower SFRs). As entropy rises with radius, the points on the leftmost panel of course also rise (i.e.. larger local entropies will be found at 
larger maximal FUV radii). Both results are found to be roughly consistent (although with large scatter) with the respective predictions by Cavagnolo et al. 
(2008) and Sharma et al. (2011), which are marked by the green dashed line on both plots. 


works to raise the local cooling time, resulting in the system’s self¬ 
regulation at nearly t coo i ~ 10 (Voit & Donahue 2015; Voit et al. 
2015a, b). 

The maximal radius within which FUV emission is detected 
in our images can be used in a very rough comparison with this 
theoretical prediction. Assuming (a) that the star formation in our 
sources is indeed powered by a precipitation-based cooling flow 
and (b) there is no non-detected star formation beyond the largest 
measured FUV radius, this “maximal FUV radius” serves as a 
rough observable tracer of the radial threshold for the onset of cool¬ 
ing flow powered star formation. Models of precipitation-regulated 
feedback predict that this radius should coincide with t coo \ & 10hr. 
Of course, assumption (b) may be not reasonable, as deeper FUV 
observations may reveal larger maximal FUV radii for a significant 
fraction of our sample. One must also consider that these filaments 
may have been uplifted by jets or buoyant cavities in some cases. 
The maximal FUV radius should therefore be treated as a lower 
limit in this caveat-laden test. 

To derive we adopt the spectrally deprojected X-ray emis- 
sivity, cooling time, and electron density profiles from work on the 
ACCEPT sample by Donahue et al. (2006) and Cavagnolo et al. 
(2009). To these, we fit third-order polynomials in log space, and 
then analytically differentiate to obtain the gravitational free fall 
time tf[. The presence of the BCG was accounted for by enforcing 
a minimum value of the gravitational acceleration g equal to that 
of an isothermal sphere with a velocity dispersion of 250 km s -1 (a 
correction that is only important at radii <^10 kpc). 

The result is shown in Fig. 13, where we compare the maximal 
FUV radius with the X-ray entropy and cooling-to-dynamical time 
ratios measured at that same cluster-centric radius. Within our sam¬ 


ple, the average maximal radius (and its ±lcr interval) within which 
FUV continuum emission is detected is 14.7 ± 10.14 kpc. At this 
average radius, the average (and ±1 a) cooling time, entropy, and 
cooling-to-dynamical time ratio is 1.09 ± 0.68 Gyr, 28.27 ± 12.12 
keV cm 2 , and 14.78 ±8.57, respectively. These average values and 
their ±lcr intervals are marked by the black solid and dashed lines 
(respectively) on the two panels in Fig. 13. Blue and red points are 
used to differentiate between those sources with unresolved radio 
sources (and higher SFRs) and resolved radio sources (and lower 
SFRs). As entropy rises with radius, the points on the leftmost panel 
of course also rise (i.e., larger local entropies will be found at larger 
maximal FUV radii). 

We find that the observed average “star formation onset” 
threshold of S = 28 keV cm 2 is very close to the empirical (rough) 
threshold of S « 30 keV cm 2 from Cavagnolo et al. (2008), which 
we mark with the green dashed line on the left panel of Fig. 13. 
The average observed tcooi/hf ratio of 14.78 is close to the pre¬ 
dicted threshold of t coo i/tf[ <; 10 by Sharma et al. (2012a), which we 
mark with the green dashed line on the rightmost panel of Fig. 13. 
Sources with higher star formation rates cluster more strongly 
around £ C ooi ~ 10than do sources with lower star formation rates. 
Heeding the strong caveats noted above, this may be roughly con¬ 
sistent (or at least not obviously inconsistent) with theoretical pre¬ 
dictions by Sharma et al. (2012a), McCourt et al. (2012), and Voit 
et al. (2015a). 

Li & Bryan (2014a,b) and Li et al. (2015) have published 
new adaptive mesh hydrodynamical simulations of BCGs consis¬ 
tent with the precipitation-regulated AGN feedback framework. In 
Fig. 14, we show selected snapshots of the cold, star forming gas 
produced by the standard run in Li et al. (2015) for a single BCG, 
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Observation Simulation Observation Simulation 


Figure 14. Observations compared with simulations by Y. Li and collaborators (e.g., Li & Bryan 2014a, b; Li et al. 2015). (Top panels ) A multiwavelength 
(X-ray, FUV, and Ha) composite of Hydra A, compared with a single snapshot of the Li et al. (2015) simulation at 2.9 Gyr. The RGB channels are set to 
roughly simulate Ha, FUV, and X-ray emmissivities, respectively, though two “cheats” have been used: the green channel shows cold, dense gas where star 
formation appears in the simulation — the assumption, then, is that FUV emission from these young stars would roughly show the same morphology. The 
red channel does not explicitly show simulated Ha, but rather intermediate temperature gas whose morphology closely matches a rough scaling to simulate 
collisionally ionized (but not photoionized) Ha. It should therefore be treated as a very rough approximation. (Bottom panels ) Selected FUV images from our 
sample are shown in orange, and projected density-weighted density snapshots from the simulation are shown in blue. The bright knots and filaments show 
the high density low temperature clouds that are forming stars in the simulation (see Li et al. 2015 for details). Our sample’s redshift range spans ~ 3 Gyr of 
cosmic history. The simulation is capable of producing star forming structures similar in axis ratio, physical extent, and star formation rate within the same 
~ 3 Gyr temporal slice of the simulated cluster’s evolution. 
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and compare these with our FUY images. The simulated cooling 
flow begins at t ~ 300 Myr (roughly the central cooling time of 
Perseus), which ignites AGN feedback. The jets trigger more ICM 
to cool into filamentary structures, causing more cold gas to rain 
into the nucleus. This fuels both star formation and black hole ac¬ 
cretion, leading to a major AGN outburst. AGN feedback heats up 
the core, reducing the cooling rate and causing Loo 1 to increase; star 
formation gradually consumes the cold gas until it vanishes, which 
turns off AGN feedback and allows the ICM to cool again. The 
cluster experiences three such cycles within 6.5 Gyr. 

Most of the images shown in Fig. 14 are from the second cy¬ 
cle that starts around t- 1.6 Gyr. Within this cycle, the simulation 
produces star forming structures that are remarkably similar to our 
FUV images not only in morphology, axis ratio, and physical ex¬ 
tent, but also in star formation rate. That is, when the simulated 
morphology matches that of one of our FUV images, its associ¬ 
ated star formation rate also roughly matches. These structures are 
reproduced by the simulation in the same ~ 3 Gyr temporal slice 
spanned by the redshift range of our sample. Their simulated mul¬ 
tiphase (hot, warm, and cold) gas morphology is also remarkably 
similar to the multiwavelength morphologies we have presented 
here. This is illustrated in the topmost panels of Fig. 14, in which 
the RGB multiwavelength composite of Hydra A in X-ray, Her, and 
FUV is nearly indistinguishable from the simulation snapshot at 2.9 
Gyr, which comes complete with a rotating, star forming disk. 

While surely other stochastic events such as mergers play an 
additional role in sculpting the morphology of star forming fila¬ 
ments in CC BCGs, the Li et al. simulations show that they need not 
be invoked to explain the FUV morphologies we observe (the same 
was true for a similar comparison made in Donahue et al. 2015). 
Instead, nearly all FUV morphologies shown in this paper appear 
in a simulated BCG whose evolution is driven by precipitation- 
regulated AGN feedback. The FUV images presented in this paper 
may be snapshots at successive stages of an ICM cooling / AGN 
heating cycle. 


5 SUMMARY & CONCLUDING REMARKS 

We have analysed the far ultraviolet morphology of star forming 
clouds and filaments in 16 low-redshift (z < 0.29) cool core bright¬ 
est cluster galaxies. X-ray, Lya, Ha, broadband optical/IR, and ra¬ 
dio maps were compared with the FUV emission, providing a high 
spatial resolution atlas of star formation locales relative to the am¬ 
bient hot and warm ionised gas phases, as well as the old stellar 
population and radio-bright AGN outflows. The main results of this 
paper are summarised as follows. 

• Nearly half of the sample possesses kpc-scale narrow fila¬ 
ments that, in projection, extend toward, into, and around radio 
lobes and/or X-ray cavities. Most (but not all) of these filaments 
are FUV-bright and forming stars, and we suggest that they have 
either been uplifted by the radio lobe or buoyant X-ray cavity, or 
have formed in situ by jet-triggered star formation or rapid cooling 
in the cavity’s compressed shell. 

• The maximal projected radius to which FUV emission is 
observed to extend corresponds to a cooling-to-freefall time of 
tcooi/tff ~ 10 for the majority of the sample. Sources with higher 
star formation rates cluster more strongly about this ratio than do 
sources with lower star formation rates. This may be roughly con¬ 
sistent (or at least not inconsistent) with theoretical predictions 
by Sharma et al. (2012a), McCourt et al. (2012), and Voit et al. 
(2015a). We nevertheless stress that maximal FUV radius is not the 


ideal tracer for the onset of purported ICM precipitation, as one 
must consider imaging depth, extinction, and other morphological 
drivers such as filament uplift. 

• The diverse range of morphology, axis ratio, spatial extent, 
and star formation rate in our FUV sample is almost entirely recov¬ 
ered in a single simulation by Li et al. (2015), demonstrating that 
galaxy-scale stochastic events such as mergers need not be invoked 
to explain the complex FUV morphologies we observe. Instead, we 
suggest that our images represent snapshots of a prototypical cool 
core BCG at many stages of its evolution. 
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Figure Al. The Abell 11 BCG (z = 0.1660). The FUV continuum image and broadband optical unsharp mask are shown in the left and right panels, respectively. 
Lyct contours are overlaid in green. The radio source is unresolved. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 4.392 x 10 -17 
ergs cm -2 A -1 electron -1 . The centroids of both panels are aligned, with east left and north up. 
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Figure A2. A multiwavelength view of the Abell 1068 BCG (z = 0.1375). X-ray, FUV continuum, and broadband optical images are shown in the top left, 
centre, and right panels, respectively. Contours of constant FUV continuum surface brightness are overlaid in blue on the X-ray panel, and in white on the 
optical panel. The white cross on the FUV panel marks the location of the unresolved radio source. The FUV colour bar can be scaled to a flux density by the 
inverse sensitivity 4.392 x 10 -17 ergs cm -2 A -1 electron -1 . The white box on the X-ray panel marks the FOV of the two rightmost panels. The bottom panel 
shows the same broadband optical image in a different colour scale, a wider FOV, and with FUV contours shown in white. The nearby companions are labeled 
“D” and “E” to correspond to the notation used by McNamara et al. (2004) in their discussion of these companions as possible gas donors. The centroids of 
all panels are aligned, with east left and north up. 



© 0000 RAS, MNRAS 000, 000-000 











24 G. R. Tremblay et al. 


electrons / sec 

0.0E + 00 2.5E-04 5.0E-04 



30 15 0 -15 -30 

Arc Seconds 

Center: R.A. 13 03 42.56 Dec -24 14 43.6 


10 5 0-5 -10 

Arc Seconds 


10 5 0 -5 -10 

Arc Seconds 


Figure A3. The Abell 1664 BCG (z = 0.1283). The X-ray, FUV continuum, and broadband optical unsharp mask images are shown in the left, centre, and 
right panels, respectively. Lya contours are overlaid in black (on the X-ray panel) and green (on the FUV and optical panels). The radio source is unresolved. 
The FUV colour bar can be scaled to a flux density by the inverse sensitivity 4.392 x 10 -17 ergs cm -2 A -1 electron -1 . The centroids of all panels are aligned, 
with east left and north up. 
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Figure A4. The Abell 1795 BCG (z = 0.0625). The X-ray, FUV continuum, and broadband optical unsharp mask images are shown in the left, centre, and 
right panels, respectively. FUV contours are overlaid in blue and green on the X-ray and optical panels, respectively. The radio source is shown in white and 
black contours on the X-ray and FUV panels, respectively. Note the remarkable spatial correspondence between the FUV continuum emission and the radio 
source. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 2.173 x 10 -17 ergs cm -2 A -1 electron -1 . The centroids of all panels are 
aligned, with east left and north up. 
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Figure A5. The Abell 1835 BCG (z = 0.2532). The X-ray unsharp, FUV continuum, and broadband optical unsharp mask images are shown in the left, centre, 
and right panels, respectively. Lya contours are overlaid in blue on the X-ray panel and in green on the FUV and optical panels. The unresolved radio source 
is overlaid in black contours on the FUV panel. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 1.360 x 10 -16 ergs cm -2 A -1 
electron -1 . The centroids of all panels are aligned, with east left and north up. 
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Figure A6. A multiwavelength view of the Abell 2199 BCG (z = 0.0302) and its remarkable radio source. An X-ray unsharp mask, FUV continuum map, and 
optical unsharp mask are shown in the top left, centre, and right panels, respectively. The double-double FR I radio source 3C 338 is shown in white contours 
on the X-ray panel, and in green contours on all other panels. Note the ~ 20 kpc-scale X-ray cavities cospatial with the radio lobes. FUV continuum contours 
are overlaid in blue and white on the X-ray and bottom-left optical panel, respectively. The FUV colour bar can be scaled to a flux density by the inverse 
sensitivity 2.713 x 10 -17 ergs cm -2 A -1 electron -1 . The black box on the X-ray panel marks the FOV of the two rightmost panels. The bottom panel shows 
wider FOVs of the broadband optical (and optical unsharp mask) images, and clearly demonstrates that the southern component of the radio source has no 
optical counterpart. As discussed in Nulsen et al. (2013), the radio source has either restarted while the host galaxy has moved north (in projection) at very 
high peculiar velocity (unlikely), or sloshing X-ray gas has pushed the relic radio source south. The centroids of all panels are aligned, with east left and north 
up. 
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Figure A7. The Abell 2597 BCG (z = 0.0821). The X-ray unsharp, FUV continuum, and broadband optical images are shown in the left, centre, and right 
panels, respectively. 330 MHz and 8.4 GHz radio contours are shown in green and white (respectively) on the X-ray panel, and 8.4 GHz contours are shown in 
black on the FUV and optical panels. Lya contours are overlaid in blue on the X-ray panel and in green on the FUV and optical panels. The FUV colour bar 
can be scaled to a flux density by the inverse sensitivity 4.392 x 10 -17 ergs cm -2 A -1 electron -1 . The multiwavelength data for A2597 are discussed at length 
in Tremblay et al. (2012a, b). The centroids of all panels are aligned, with east left and north up. 
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Figure A8. NGC 4696, the BCG of the Centaurus Cluster (z = 0.0099). The X-ray surface brightness map is shown in the leftmost panel, with radio and FUV 
contours overlaid in white and blue contours, respectively. The centre-top panel shows the 5-band optical map with radio contours overlaid in green. The 
centre-bottom panel shows a B/H-bmd colour map that highlights the dramatic (and well-known) 5 kpc-scale dust lane. FUV contours are overlaid in blue. 
The rightmost panel shows the FUV continuum data, with radio contours overlaid in green. The FUV colour bar can be scaled to a flux density by the inverse 
sensitivity 4.392 x 10 -17 ergs cm -2 A -1 electron -1 . The centroids of all panels are aligned, with east left and north up. 
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Figure A9. A multiwavelength view of the Hydra A (Abell 780) BCG (z = 0.0549), one of the most mechanically powerful radio sources in the known Universe 
(Wise et al. 2007). An X-ray unsharp mask, FUV continuum map, and optical unsharp mask are shown in the left, centre, and right panels, respectively. The 
FRI radio source 3C 218 is shown in green contours on all panels. Note the ~ 50 kpc-scale X-ray cavities cospatial with the radio lobes. The FUV colour bar 
can be scaled to a flux density by the inverse sensitivity 2.713 x 10 -17 ergs cm -2 A -1 electron -1 . The white box on the X-ray panel marks the FOV of the two 
rightmost panels. The centroids of all panels are aligned, with east left and north up. 
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Figure A10. NGC 1275, the BCG of the Perseus Cluster (z = 0.0176). The X-ray, Ha+[N II], and FUV continuum maps are shown in the left, centre, and right 
panels, respectively. Radio contours are shown in green on the X-ray panel, and FUV contours are shown in blue and black on the X-ray and optical panels, 
respectively. The white box on each panel shows the FOV of the panel to the right. Some of the FUV emission may be attributable to an unrelated galaxy that 
is superimposed on the line of sight. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 2.713 x 10 -17 ergs cm -2 A -1 electron -1 . The 
centroids of all panels are aligned, with east left and north up. 
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Figure All. A multiwavelength view of the PKS 0745-191 BCG (z = 0.1028). The two leftmost panels show an X-ray surface brightness map at top and a 
X-ray unsharp mask at bottom (whose zoomed-in FOV is marked by the white box on the top panel). Radio and FUV contours are overlaid in black and blue 
on the top panel, and in black and white on the bottom panel, respectively. The center and rightmost panel shows the FUV continuum map and an optical 
unsharp mask image, with radio and FUV contours overlaid in green and blue, respectively. The FUV colour bar can be scaled to a flux density by the inverse 
sensitivity 2.713 x 10 -17 ergs cm -2 A -1 electron -1 . Note that the FUV continuum follows the “spine” of the “bird-shaped” radio source. The white box on the 
bottom-most X-ray panel marks the FOV of the two rightmost panels. The centroids of all panels are aligned, with east left and north up. 
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Figure A12. A multiwavelength view of the RX J1504.1-0248 BCG (z = 0.2153). X-ray, FUV continuum, and optical unsharp mask images are shown in the 
the left, centre, and right panels, respectively. FUV continuum contours are shown in blue on the X-ray panel, and the unresolved radio source is shown in 
white contours on the FUV panel. Note the ~ 20 kpc stellar filament (FUV and broadband optical) extending along the BCG major axis. The FUV colour bar 
can be scaled to a flux density by the inverse sensitivity 1.360 x 10 -16 ergs cm -2 A -1 electron -1 . The white box on the X-ray panel marks the FOV of the two 
rightmost panels. The centroids of all panels are aligned, with east left and north up. 
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Figure A13. The RX J2129.6+0005 BCG (z = 0.235). The X-ray, FUV continuum, and broadband optical images are shown in the left, centre, and right panels, 
respectively. Lya contours are overlaid in black on the X-ray panel, and in green on the FUV and optical panels. The radio source is unresolved. The white box 
on the X-ray panel marks the FOV of the FUV and optical panels. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 1.360 x 10 -16 
ergs cm -2 A -1 electron -1 . The centroids of all panels are aligned, with east left and north up. 
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Figure A14. The ZwCl 0348 BCG (z = 0.255). The X-ray, FUV continuum, and broadband optical images are shown in the left, centre, and right panels, 
respectively. Lyo contours are overlaid in white on the X-ray panel, and in green on the FUV and optical panels. The radio source is unresolved. The white box 
on the X-ray panel marks the FOV of the FUV and optical panels. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 1.360 x 10 -16 
ergs cm -2 A -1 electron -1 . The centroids of all panels are aligned, with east left and north up. 
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Figure A15. The ZwCl 3146 BCG (z = 0.291). The X-ray, FUV continuum, and broadband optical images are shown in the left, centre, and right panels, 
respectively. Lya contours are overlaid in blue on the X-ray panel, in green on the FUV panel, and in blue on the optical panel. We also overlay FUV contours 
are in green on the optical panel. The radio source is unresolved. The white box on the X-ray panel marks the FOV of the FUV and optical panels. The FUV 
colour bar can be scaled to a flux density by the inverse sensitivity 1.360 x 10 -16 ergs cm -2 A -1 electron -1 . The centroids of all panels are aligned, with east 
left and north up. 
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Figure A16. The Zw8193 BCG (z = 0.1829) may be undergoing a (minor?) merger. FUV continuum and a broadband optical unsharp mask are shown in 
the left and right panels, respectively. Lya contours are overlaid in green (innermost contours have been removed to aid viewing). Note the highly disturbed 
morphology of the optical counterpart, as well as its apparent double nucleus. Brighter FUV emission associated with ongoing star formation is cospatial 
with the northernmost optical nucleus. Clumpy tendrils wind counter-clockwise from this bright northern FUV knot, suggestive of non-negligible net angular 
momentum perhaps stemming from the merger. The FUV colour bar can be scaled to a flux density by the inverse sensitivity 4.392 x 10 -17 ergs cm -2 A -1 
electron -1 . The centroids of both panels are aligned, with east left and north up. The radio source in Zw8193 is unresolved. 
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